Quantum effects in atomically perfect specular spin valve structures 
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A simple tight-binding theoretical model is proposed for spin dependent, current-in-plane trans- 
port in highly coherent spin valve structures under specularity conditions. Using quantum- 
mechanically coherent and spatially quantized Fermi states in the considered multilayered system, a 
system of partial Boltzmann kinetic equations is built for relevant subbands to yield the expressions 
for conductance in parallel or antiparallel spin valve states and thus for the magneto-conductance. 
It is shown that specularity favors the magnetoresistance to reach its theoretical maximum for this 
structure close to 100%. This result is practically independent of the model parameters, in particular 
it does not even need that lifetimes of majority and minority carriers be different (as necessary for 
the quasiclassical regimes). The main MR effect in the considered limit is due to the transformation 
of coherent quantum states, induced by the relative rotation of magnetization in the FM layers. 
Numerical calculation based on the specific Boltzmann equation with an account of spin-dependent 
specular reflection at the interfaces is also performed for a typical choice of material parameters. 

PACS numbers: 75.10.Hk; 75.30. Gw: 75.70.Cn; 76.50.+K 
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I. INTRODUCTION 

Fabrication of novel nanostructured spintronics devices 
and the related experimental studies of spin-dependent 
electron transport stimulate new theoretical approaches 
to the physical properties of nanosystems where quan- 
tum coherence effects can have a decisive role, in contrast 
to the mostly quasiclassical framework of the traditional 
electronics. One important class of such systems con- 
cerns spin valves [l[ formed by two ferromagnetic (FM) 
layers separated by a thin non-magnetic (NM) spacer. 
The magnetization of one of the FM layers (called pinned 
layer) is fixed by the bias from underlying antiferromag- 
netic (AFM) layer, while the magnetization of the other 
FM layer (free layer) easily rotates when a small mag- 
netic field is applied. This significantly affects the in- 
plane conductance, leading to relatively high MR values, 
typical for giant magnetoresistance (GMR) [§], but the 
technology still demands further improvements. One of 
them consists in the introduction of nano-oxide layers 
(NOL's) just above the free layer and inside the pinned 
layer (so that the pinning is not disrupted) Such 
NOL-equipped device, the so called specular spin valve 
(SSV, Fig. [TJd) can more than double the GMR ratio of 
simpler stacks (Fig. Hk). The increase of MR is believed 
to arise from the specular reflection of electrons at the 
FM/NOL interfaces. 

But, besides the evident effect of carriers confinement, 
the reduced normal-to-plane scale d of magnetic layers 
(few nm thickness, controlled within a 1 A precision) 
might allow for a pronounced quantization of the nor- 
mal component of quasi-momentum, as already indicated 
by the recent data on spin-resolved electronic reflection 
from magnetic nanolayers Q, p]- Furthermore, it is ex- 
pected that the relevant modes at the Fermi level for each 
polarization are dramatically restructured when the mu- 
tual polarization of magnetic layers is changed. All this 
can qualitatively change the kinetics of spin-dependent 
transport, compared to the usual diffusive scenario for a 
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Figure 1: Schematics of spin valve structures: a) common and 
b) specular. 



quasi-continuous spectrum However, the microscopic 
understanding of electron specular interface reflection is 
still far from complete, in particular its role in size quan- 
tization and coherence of Fermi states. Here we propose 
a theoretical description of these effects, through a prop- 
erly modified Boltzmann kinetic equation, taking into 
account the formation of transverse-quantized electronic 
subbands and spin-dependent specular reflection at the 
interfaces within the simplest tight-binding model, easy 
enough to advance up to numerical calculations of the 
MR behavior. 



II. MODEL 

Let us begin from a single metal layer, made of n 
atomic planes with simple cubic lattice coordination and 
hopping integral t between nearest neighbors at distance 
a. The respective electronic spectrum for given spin po- 
larization a =|, I and planar quasimomentum k consists 
of n subbands of the form £ a ,k,<r = £k + A CT + 5 a (Fig. 
[21). Here = 2t (2 — cosak x — cosak y ) is the 2D disper- 
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Figure 2: Sketch of the dispersion laws (along the diagonal 
k x = k y = kj_ of 2D Brillouin zone) in spin-splitted and spa- 
tially quantized subbands of a magnetic nanolayer. The cir- 
cles indicate Fermi momenta for particular (minority spin) 
subbands and the related Fermi velocities v a correspond to 
the slopes of dispersion laws. 



sion law for a single plane, and in a ferromagnetic metal 
it is accompanied by the Stoner energy shift A CT = ±A 
for minority and majority spins respectively. The spa- 
tial quantization is accounted for by the subband shifts 
6 a (a — 1, ... ,n) which are the eigenvalues of the n x n 
secular equation 
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with exact values S a = 2tcos[ira/(n + I)]. The wave 
function for the a, k, a state, at the planar position r 
in the j-th plane, is ^> a ,k,a-(r,i) = Aj- e lk ' r x CT , where 
the components of the n-dimensional eigenvector A^ a ' 
related to the eigenvalue S a are explicitly given by 
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and Xct is the spin function. 

Next the model is extended to include the hopping 
t' between the neighbor FM and NM layers, hybridiz- 
ing the subbands of the free FM layer (composed of nf 



atomic planes) 
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of the NM spacer (composed of 
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Figure 3: Energy band structure in the trilayered system. 
All the modes are doubly degenerate and the relevant modes 
at the Fermi level are marked with circles. Inset: spatial 
composition of atomic planes forming the sets Jf, s , P in /-, s-, 
and p-layers. 



n s planes) e s a k (with A = 0), and of the pinned FM 
layer (n p planes) £„ kcr . We shall denote the respec- 
tive eigenvectors (for the uncoupled layers) by F^ a \ S^ a \ 
and P( a \ with the components given again by Eq. [2] 
for n — nj,n p ,n s , while the notation is adopted 

for the eigenvectors of the coupled system. The spec- 
ularity effect in this approach is modeled by zero cou- 
pling of the FM layers to their outer neighbors. The 
resulting spectrum totals up to nt = 2(n/ + n s + n p ) 
spin-resolved modes with energies e Q ,k and wave func- 
tions *a,k(r,j) = M^ a) e lk ' r Xa( a ), where a = l,...,n t 
and a(a) is the implicit polarization of a-th mode (Fig. 
[3]). We emphasize that from the total of n t modes, only 



a smaller number, n r , of modes, those present on the 
Fermi level, are relevant for conductance. Thus, for the 
characteristic case of FM Co layers, only minority spin 
subbands should take part in the transport (as suggested 
by the bulk Co band structure j§|). Moreover, we have 
to take into account the sizeable differences in the corre- 
sponding Fermi velocities v a (practically coincident with 
those in the uncoupled layers, Fig. [2]). The most es- 
sential effect of hybridization is on the amplitudes 
which are generally some weighted combinations of all the 
F, P, S modes, and the crucial point is that the weights 
of F, P components in the relevant modes are strongly 
dependent on the mutual polarization of FM layers (see 
below) . 

Then the kinetics of the composite system is described 
by the set of n r distribution functions / aj k = /^° k + g aj k 

where /^° k = [e /3 ( £a k_£F ) + l] 1 is the usual equilibrium 
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J (ef — £/3,k), that is to the average of the non- 

equilibrium distribution over the Fermi surface and so it 
should vanish. Then the solution takes the common form 
9a, k = h^ 1 r a e'E ■ df^/dk where the relaxation time for 
the a-th mode is defined by 
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including the Fermi densitiy of states pp = 
(a/2ir) 2 f dkS — ep) for each (3-th mode. Then 
the total conductivity is found from Eq. [3] as a sum of 
partial contributions: 
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Figure 4: Configurations of Fermi lines for spatially quantized 
subbands of minority electrons in the Brillouin zone. The 
characteristic points along high symmetry directions YM and 
MX were used to approximate the averages of n^ 1 and v 2 a . 



Fermi function with j3 = \/k^T and g a ^ is the non- 
equilibrium part due to the external electric field E. The 
current density is given by the sum 
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where means summation over the n r relevant modes, 
v Q ,k = ?i _1 <9£ Q ,k/<9k is the electron velocity, and the 
components of the non-equilibrium distribution are de- 
fined from the system of Boltzmann equations: 
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?a,k) = 0. (4) 



Here J^" means summation over relevant modes with 

conserved spin, cr(a) — a(f3), and is the transition 

rate due to scattering from k state of a-th subband to 
k' state of /3-th subband. We consider transitions only 
due to elastic scattering by random point-like impurities 
with potential V and concentration c <C 1 (per unit cell) . 
Then the Fermi Golden Rule transition rates are ui^'k = 
^a,/3<5 (Sa.k — £/3,k') with the scattering factors (averaged 
in impurity positions) 
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In this simple model, the first term in the collision 
integral of Eq. [4] turns out to be proportional to 



where (v%) ~ Pa 1 (a/^n) 2 J dkv^ k S (e Qjk - e F ) is the 
average of the respective squared Fermi velocity. In fact, 
this is a particular case of the general Landauer formula 
[13], written for the present system of n r coherent quan- 
tum channels. 

The system, Eqs. \VjjE[ can be routinely treated by nu- 
merical methods at any relative orientation of magnetiza- 
tions in /- and p- layers, from parallel (TT) to antiparallel 
(T!), to result in the principal quantity of interest, the 
magnetoresistance 
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But some qualitative conclusions about the specularity 
effect on MR in a nanolayered device can be drawn al- 
ready from simple inspection of the discrete structure of 
the amplitudes M- , according to the following remarks. 

First of all, we suppose that in absence of hybridization 
the majority and minority subbands are well separated 
from each other and from the spacer subbands (like the 
situation in bulk Co and Cu). Then we notice that the 
j-configurations of the above amplitudes are essentially 
different for |T and || cases and hence consider them 
separately. Finally, an important factor for the very ex- 
istence of GMR (in this quantum coherent conductance 
regime) is the presence of certain "resonances" between 
relevant modes at the Fermi level. Namely, a resonance 
appears between two (unhybridized) modes £^ k a and 
e /3 k a ^ their energy separation near the Fermi level is 
less then the effective coupling ~ t' 2 je s (mediated by the 
spacer modes at typical energy distance e s , see Fig. [3]). 
Moreover, for the sake of clarity, we shall restrict the fol- 
lowing consideration to the simplest situation of identical 
/- and p-layers where all rif — n p modes are relevant and 
can form resonant /p-pairs. 

Thus, in the TT configuration, there appears a strong 
hybridization in each F a , P a pair, forming two collective 
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modes as their bonding and anti-bonding combinations 
(in neglect of small contributions ~ t' 2 / (e s A) <C 1 of the 
rest of the modes): 
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where Jf, s ,p are the sets of atomic planes entering /-, 
s-, and p-layers (see inset in Fig. [3]). The respective 
relaxation times are given by 
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Then we can use the exact sum rule for the amplitudes, 
Eq. H 



E(^) 2 = ^r(i + ^ + 2- n+1 " g )' 
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to present the relaxation times, Eq. \W[ as 
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Contrariwise, in the || configuration, all the relevant 
modes remain almost unhybridized, taking nearly "local" 
forms: 
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and in this approximation we obtain for the relaxation 
times \ half the value of Eq. [121 Then the magnetore- 
sistance, Eq. [8l is readily estimated as AR/R w 100%. 
We notice that this result is practically independent of 
the parameters of interlayer coupling and impurity scat- 
tering, in particular it does not even need that lifetimes 
of majority and minority carriers be different (as nec- 
essary for quasiclassical regimes). The main MR effect 
in the considered limit is due to the variation of coher- 
ent quantum states, induced by the relative rotation of 
magnetization of the FM layers. 



III. NUMERICAL CALCULATIONS 

To certify the above qualitative considerations, a de- 
tailed numerical calculation was done for the particular 
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= i TT 


0.5175 


0.3431 


0.3417 


0.5151 
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1.0347 


0.6857 


0.6824 


1.0287 
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= 2 TT 


0.3431 


0.5118 


0.5097 


0.3415 
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■i 


0.6857 


1.0224 


1.0175 


0.6818 
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= 3 TT 


0.3417 


0.5097 


0.5077 


0.3402 




Tl 


0.6824 


1.0175 


1.0127 


0.6786 


a 


= 4 TT 


0.5151 


0.3415 


0.3402 


0.5127 




Tl 


1.0287 


0.6818 


0.6786 


1.0228 



Table I: Scattering factors (Cl a ,@) for the Fermi modes (a, (3 - 
1, 2, 3 and 4) in the parallel (TT) and antiparallel (Tl) config 
urations. 
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10 w p a {3- 1 ) 


0.209 


0.2439 


0.3238 


0.5517 


10 10 (^>(m 2 /s 2 ) 


0.7431 


2.2595 


2.9334 


2.6008 


10- 12 rP( s) 


1.7018 


1.8157 


1.823 


1.7094 


1(T 12 tP(s) 


0.8515 


0.9091 


0.9133 


0.8563 



Table II: The Fermi density of states p a , averages of squared 
Fermi velocities (w 2 ) and relaxation times r a for the Fermi 
modes (a = 1, 2, 3 and 4) in the parallel (TT) an d antiparallel 
(Tl) configurations. 



choice of parameters: t = t' = 0.25 eV, A = 0.5 eV, 
e s = 2 eV (a single-band model for real d-bands of Co 
and Cu), nf = n p = 4, n s = 3 (a simple discrete struc- 
ture of layers), and V = 0.5 eV, and c = 0.01 (typical 
impurity parameters). The Fermi velocities (and their 
inverse values) for two characteristic directions in Bril- 
louin zone were used to approximate the partial densities 
of states 
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(L a being the length of respective Fermi line, Fig. [4]), 
and then (v a ) w v a ,rAiv a ,MX ■ The obtained numerical 
results for fl a ,(3, p a , T a , and (v a ) are illustrated in tables 
U and [ill respectively, for the TT and Tl configurations. 
These numerical values lead to w 99.65% of magnetore- 
sistance, that is quite close to the maximum possible MR 
= 100% in the coherent regime. To compare with, for 
purely incoherent currents there will be no MR at all 
in such two-layer system, so that the finite effect only 
appears from their partial mixing due to scattering at 
the interfaces and is estimated phenomenologically as 
~ exp (—d/l) of the above maximum value. 

Actually, the experimental MR values in specular spin 
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valves [lll.[l2| are clearly lower than the above model esti- 
mates. This can be due to a number of important factors, 
not included into the present simple model (which there- 
fore should be considered as a certain idealized reference 
case). First of all, the postulated ideal specularity condi- 
tion (supposing the wave function fully confined within 
the n-plane system) cannot be exact in reality, and a 
considerable part of the electronic density can "escape" 
through the NOL barriers to adjacent non-magnetic (or 
AFM) layers. This part would act as a parallel con- 
duction channel, mostly unchanged at reorientation of 
FM electrodes and hence restricting the magnetoresis- 
tive effect. Also, the used model of rigid Stoner shifts of 
spin subbands in FM electrodes of course overestimates 
sharpness of the spin-dependent energy barrier between 
these electrodes, where in fact the band structures are 
not uniform on scales of few atomic layers. Other re- 
strictive factors are the temperature effect (by phonons 
and magnons), the roughness in the FM/NM interfaces 
and the presence of defects as grain boundaries, dis- 
placements and distortions in the crystalline structure, 
which all reduce the coherence of relevant quantum states 
and so the validity of the Landauer formula. At last, 
the single-band model may be oversimplified, compared 
to the real hybridized s-d band structures of transition 
metals used in various numerical studies of spin-valves 
[13, EH, EH, but it would be rather problematic to keep 
the analytic description at such detailed level. Neverthe- 
less, a further development within the present model can 



be done, varying the number n r of relevant modes and 
admitting the presence of both spin polarizations among 
these modes. 



IV. CONCLUSIONS 

A simple single-band tight-binding model was devel- 
oped to estimate theoretically the maximum possible en- 
hancement of GMR in the system of quantum coherent 
FM nanolayers, using a specific set of Boltzmann equa- 
tions for spatially quantized and spin-resolved subbands 
and a Landauer-type formula for the spin-dependent con- 
ductance. It is shown that a limit GMR value close to 
100% can be reached for a fully coherent (and fully spec- 
ular) SSV nanostructure and the reducing factors for this 
value in real SSV systems are discussed. 
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